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Abstract. The relation between the properties of a specific crystallographic site and 
the properties of the full crystal is discussed by using spherical tensors. The concept 
of spherical tensors is introduced and the way it transforms under the symmetry 
operations of the site and from site to site is described in detail. The law of spherical 
tensor coupling is given and illustrated with the example of the electric dipole and 
quadrupole transitions in x-ray absorption spectroscopy. The main application of the 
formalism is the reduction of computation time in the calculation of the properties 
of crystals by band structure methods. The general approach is illustrated by the 
examples of substitutional chromium in spinel and substitutional vanadium in garnet. 

Uncomment for PACS numbers title message 



PACS numbers: 61.72.S-,61.05.cj,02.20.-a,61.50.Ah 



Site symmetry and crystal symmetry 



2 



1. Introduction 

This paper deals with the relation between the property of specific sites in a crystal 
and the property of the crystal itself. We shall discuss in particular the case of x-ray 
absorption spectroscopy but many of our results are general. 

The aim of this paper is to introduce x-ray physicists, experimentalists and 
theoreticians, to the use of spherical tensors to analyze theoretical and experimental 
spectra. Therefore, our presentation is as pedagogical as we can manage. The 
manipulation of spherical tensors can be quite complicated, we have tried to reduce 
it to the simplest possible rules. 

An atom in a crystal occupies generally a crystallographic site that is not unique. 
The present paper was inspired by the concrete case of a chromium atom sustituted 
for aluminum in spinel pp. The chromium atom can occupy sixteen equivalent trigonal 
sites in the cubic cell. The (normalized) x-ray absorption spectrum of chromium is 
the average of the spectra of the chromium atoms occupying the sixteen sites. The 
angular-dependent spectrum of a chromium atom in a specific site is determined by 
its environment, but some features of this spectrum are lost by the average over the 
equivalent sites. If we want to investigate the deformation of the environment due to 
the substitution of chromium for aluminium, we need to know the spectrum of chromium 
at a given site. But some features of this spectrum are washed out by the fact that we can 
only measure an average spectrum. The present paper describes precisely the relation 
between the spectrum of a given site and the spectrum of the crystal. 

If we want now to calculate the spectrum of chromium, we need to put a chrominum 
atom at a specific aluminum site, relax the environment, carry out a self-consistent 
calculation with and without core hole, and calculate the spectrum of this site. Then, 
we need to repeat the procedure for all equivalent sites and take the average spectrum. 
A considerable amount of time can be saved if we calculate the spectrum of one site and 
deduce the spectrum of the other sites by symmetry considerations. The present paper 
gives all the tools required to do so. 

Let us first defend the case of spherical tensors, that will be precisely defined in 
the next section. Most physical properties are tensors and they are usually taken to be 
Cartesian tensors. In this paper, we want to stress the advantage of using spherical 
tensors. They are a refinement with respect to Cartesian tensor. For example, a 
second-rank Cartesian tensor is the sum a zeroth-rank, a first-rank and a second-rank 
spherical tensors. This refinement enables us to discard irrelevant contributions. For 
example, the quadrupole transitions in x-ray absorption spectroscopy are represented 
by a fourth-rank Cartesian tensor with 81 components, whereas it is represented by 
a zeroth-rank, a second-rank and a fourth-rank spherical tensor with 15 components 
altogether. Moreover, the spherical average used to represent powder sample is obtained 
by taking only one spherical tensor component (the zeroth-rank tensor) whereas it is the 
sum of twelve Cartesian components. Thus, using spherical tensors can save computing 
time. 
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The drawback of the spherical tensor analysis is that some subtelties enter its 
practical use (the same is true for Cartesian tensors). A purpose of this paper is to give 
a detailed presentation of these problems and of their solutions. 

In practice, one faces the frustrating task of sorting out the various conventions 
used by the authors and to determine whether the rotation is active or passive, if the 
row index of the Wigner matrix increases from left to right as in a normal matrix or 
decrease as in reference [2], p. 47, which definition of the Wigner matrix is used (seven 
inequivalent definitions are found in the literature [3]), how the tensor transforms under 
complex conjugation, etc. Misprints are also numerous in the literature and we made 
thorough efforts to give correct formulas. 

The paper starts with a presentation of the Cartesian and spherical tensors and a 
detailed description of the notation. Then, the concept of tensor coupling is introduced 
and simple formulas are given to change the coupling order. These methods are applied 
to the calculation of electric dipole and quadrupole transitions. This completes the 
generalities on spherical tensors. Then, we consider the case of crystals. We show 
how the symmetries of a crystal site constrains the spherical tensor components in the 
reference frame of the site and in the reference frame of the crystal. We also give the 
precise formula relating these reference frames. As an illustration, we treat the example 
of the Al site in spinel and garnet. Finally, we determine the spherical tensor describing 
the full crystal from those of the sites. An appendix gathers the formulas that were 
used to calculate rotation matrices, Wigner matrices and solid harmonics. 

2. Cartesian and spherical tensors 

A Cartesian tensor is the generalization of a scalar, a vector or a matrix. A vector r = 
(x, y, z) is transformed by a rotation R into a vector r' = Rr, so that = Xy=i -^*i r i- 
A vector is also called a first-rank tensor. A matrix A defined by Ay = r^r,- transforms 
under rotation into 



Any matrix Ay that transforms under rotation as in equation (TjQ) is called a second-rank 
Cartesian tensor. More generally, an nth-rank Cartesian tensor is a generalized matrix 
Ai 1 _i n that transforms under rotation as in equation ([I]), but with the product of n 
matrices R instead of just two. Cartesian tensors are ubiquitous in physics. As we shall 
see, the electric dipole transitions are described by a second-rank tensor, the quadrupole 
ones by a fourth-rank tensor. 

Cartesian tensors have very simple transformation rules under rotation, but they 
suffer from a severe drawback: they are not irreducible. To see what that means, consider 
a second-rank tensor Ay. Its trace is t — X^=i A-u. It transforms under rotation into 




(1) 



Im 
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where we used the fact that RR T = Id, where R T is the transpose of R. We recover 
the fact that the trace of a matrix is invariant: it transforms into itself by rotation. A 
second-order Cartesian tensor gives another interesting object, the vector v, which is 
defined by = Yljk e ijkAjk, where the sum is over j and k from 1 to 3. The Levi- 
Civita symbol is 1 if (i,j,k) is a cyclic permutation of (1,2,3), it is -1 if (i,j,k) 
is another permutation of (1,2,3) and it is if two indices are identical. Then, using 
the identity J2jk e ijkRjiRkm = Ylj Rij e jim, we see that v transforms under rotation as 
a vector: v' = -Rv. Therefore, from a second-rank Cartesian tensor, we can build a 
linear combination of its elements (the trace) that is invariant under rotation, and three 
linear combinations of its elements that transform into each other as the components 
of a vector. More generally, a tensor is said to be reducible when there are linear 
combinations of its elements that transform into each other under rotation. When a 
tensor is not reducible, it is called irreducible. Thus, a vector is irreducible but a second- 
rank Cartesian tensor is reducible. The irreducible tensors are called spherical tensors 
and will be the main topic of this paper. 

The first spherical tensors were the spherical harmonics Y™. For each £, there 
are 2£ + 1 spherical harmonics Yf 1 that transform into each other by rotation. More 
precisely, for each rotation R, there is a unitary matrix D e , called a Wigner matrix (to 
be precisely defined in the next section), such that the rotation of Y™ by R is 

£ 

RY™ = Yl YfDi, m {R). 

m'=-£ 

Spherical tensors are defined in analogy with spherical harmonics. An £th rank 
spherical tensor, denoted by T e , is a set of 2£ + 1 components, written where 
m = —£, —£ + 1, ...,£ — 1,£, that transform under rotation as 

t 

RTfn = T^/D^^R). (2) 

m'=—£ 

This definition is not very concrete, but we shall see how spherical tensors are built in 
practice. It is an unfortunate but historical fact that the position of £ and m is different 
in the spherical harmonics Y™ and the spherical tensors T^. 

When many spherical tensors are involved in a formula, we use also the notation P a , 
Q b , etc. For notational convenience, we shall often write RT e = T e D e (R) for equation 
(J2J), as for the product of a matrix and a vector. Moreover, the product of two Wigner 
matrices will be written D e (R)D e (R'). 

2.1. Further symmetries 

Spherical tensors can satisfy other symmetries. For example, most of them are built 
from Hermitian operators and satisfy time-reversal symmetry: (see equation (4) p. 61 

of m- 

(T e m T = (-l) m Tf m . (3) 
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In this paper, the only tensors that do not satisfy time-reversal symmetry are those 
built from the polarization vector e that can possibly be complex. In that case we have 
(T^(e))* = (— l) m Tf m (e*). This happens for instance when T l is a solid harmonics (see 
section 18.21) . 

Inversion is another common symmetry. In pure electric dipole and quadrupole 
transitions inversion does not play any role because all tensors are even (i.e. invariant 
under inversion). However, the interference between dipole and quadrupole transitions 
is always an odd spherical tensor [U EJ E] . 



3. Notation 



The many conventions and misprints found in the literature forces us to describe in 
detail our notation and conventions. 

We consider active rotation, that is rotation that move the points and not the 
reference frame. For example the rotation through an angle if) about the z-axis is 
represented by 

' cos if) — sin if> ' 
R z {.ip) = sin-0 cos if) 

V l ) 

After an active rotation R, the coordinates (ri,r 2 , r 3 ) of the vector r are transformed 
into the coordinates = ^ ■ RijTj of r' = Rr. In a passive rotation, the reference frame 
is rotated: the basis vectors are transformed into the basis vectors = Y2jRij e j- 
Thus, the coordinates of a point r are transformed by the inverse matrix: r' = R~ 1 r. 

To describe the transformation of the properties of a crystal under rotation, we 
consider the case of its charge density p(r). After a rotation changing r into r' = Rr, 
the charge density p is transformed into a "rotated" charge density p' of the rotated 
crystal. To determine p', we require that the value of the charge density is not modified 
by the rotation. More precisely, we want that p'(r') = p(r). Therefore, the rotated 
function p' is defined by p'(r') = p(R~ 1 r / ). For later convenience, we denote the rotated 
function p' by Rp. The use of the same symbol R to denote the rotation of both the 
vectors and the functions should not bring too much confusion. The presence of the 
inverse rotation in the definition of Rp ensures that R'(Rp) = (R'R)p (see Ref. [7] 
p. 59). 

3.1. Wigner rotation matrices 

We denote by D l m , m (R) the Wigner rotation matrix corresponding to the rotation R. 
For example, D l m , m {R z {ifj)) = 8 mm ieT im ^ . The Wigner rotation matrices define a unitary 
representation of the rotation group, so that 

Dl^R- 1 ) = (D e mm ,(R)T, (4) 
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and 

e 

m"=-t 

3.2. Spherical harmonics 

The spherical harmonics are defined by (see Ref. [2] p. 68) 



IP + 1 

yr(OA) = ^^ r (DURe,)T, (6) 

where Rg^ is the rotation defined by the Euler angles (<ft, 6,0). For notational 
convenience, we denote by n the vector (sin 9 cos 0, sin 9 sin 0, cos 9) and we write Y^ m (n) 
and R n for Y™(0, 4>) and Rq$. This notation is justified by the fact that Y^ija) can be 
defined for any vector r, not necessarily normalized. The resulting functions are called 
solid harmonics and are described in section I8.21 Solid harmonics are required, for 
example, in the presence of elliptically polarized x-rays because n has then complex 
coordinates. A three-dimensional Cartesian vector r = (x, y, z) can be turned into a set 
of three solid harmonics 

Yr 1 (r) = J±(x-iy), (7) 

Y}(r) = -JI-(x + iy). (9) 
This definition of spherical harmonics implies 
Yr(Rn)=J2Y i m '(n)D e m/m (R- 1 ). 

m' 

This is proved by noticing that the argument Rn of the spherical harmonics corresponds 
to the argument RR n of the Wigner matrix in equation fl§]). By equation (TJJ, we 
have (D e m0 (RR n ))* = D l Qm ((RR n )~ l ) . The result follows from {RR n )~ l = R^R' 1 and 
the group representation property, equation (jSJ). The same property is true for solid 
harmonics. Therefore, 

(iar)(r) = YriR-'r) = J2 Y r'(r)Dl> m (R)- (10) 

fit' 

The presence of the spherical harmonics on the left of the Wigner rotation matrices 
ensures that R'(RY e m ) = (R'R)Y™. To show this, multiply equation (02} on the left by 
R': 

(R'(RY £ ))(r) = {R%){r)D\R) = Y t {v)D%R!)D\R) 
= Y,{v)D\R'R) = (R'R)Yi(r). 

In the foregoing proof, we simplified the notation by omitting the component index m 
(as described at the end of section [2]). 
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4. Building tensor operators 

Physical properties can be represented by spherical tensors that can often be built by 
coupling lower rank tensors. We illustrate this construction by the example of electric 
dipole and quadrupole transitions. We shall use the remarkable toolbox for spherical 
tensor calculations elaborated by Varshalovich, Moskalev and Khersonskii [3]. 

Many spherical tensors used in physics are obtained by coupling vectors. A three- 
dimensional Cartesian vector v = (x, y, z) can be turned into a first-rank spherical 
tensor v 1 by defining 

vii = (x- iy)/V2, 

v o = 

v{ = -{x + iy)/V2. (11) 

Note that solid harmonics Y £ m (v) are also spherical tensors built from v and that 
v 1 = Yi(v)-\/47r/3. However, the factor a/47t/3 is cumbersome and the definition 
v £ = Yt(v) ^/Att] {21 + 1) is often preferred. 

Two spherical tensors P a and Q b of ranks a and b can be coupled into a spherical 
tensor of rank c, denoted by {P a ^Q b } C i and defined by 

a b 

{p a ®Q%= Y,^ aah ^ pa A 

a=—a f3=—b 

where (aab/3\ay) are Clebsch-Gordan coefficients [21 [3]. The Clebsch-Gordan coefficient 
is zero if 7 7^ a + f3. The triangle relation implies \a — b\ < c < a + b. The notation 
{P a ®Q b } c is inspired by ref. [3], except for the fact that they write the rank as an 
index instead of an exponent. For example, the coupling of two vectors (i.e. first-rank 
spherical tensors) gives a zeroth-rank, a first-rank and a second-rank spherical tensors. 
The zeroth-rank tensor obtained by coupling two vectors is proportional to the scalar 
product of these vectors: {u 1 ®v 1 } = — u ■ v/y3, because (lad — «|00) = — (— l) a /\/3- 
More generally, we define the scalar product of two spherical tensors P a and Q a of the 
same rank to be (see [3], p. 64 and 65) 

a 

P a -Q a = Yl (-^P-aQa = (-l) a V2a7+T{P%Q a }°. (12) 

a=— a 

It is often necessary to modify the coupling order of the tensors. For example, 
when calculating electric dipole transitions, we have to calculate |(/|e • r|i)| 2 , where r 
is coupled to e by the scalar product and then multiplied by its complex conjugate. As 
we shall see in the next section, it is more convenient from the physical point of view 
to couple the x-ray polarization vectors e and e*. For that purpose, there is a powerful 
identity 

J pa n a X rnd Q d X V"V 1 Ng j/P°gj^ - {Q a ®S d } 9 

{P «Q } ■ {R. «S } = ^(-1) V(2a+1)i2d+1) ■ (13) 
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where g runs from \a — d\ to a + d by the triangle relation. This identity is proved in 
section O 

In the next two sections, we illustrate the recoupling methods with the calculation 
of electric dipole and quadrupole transitions. Similar methods were used to investigate 
the interference of electric and quadrupole transitions and their interferences (5J El IE] or 
to calculate x-ray scattering cross-sections [9]. 

4-1. Dipole 

The electric dipole transitions amplitudes are given by the formula Tfi = (f\e ■ r\i). 
If we denote (/|r|i) by Tfi, equation (fT2|) gives us Tfi = e ■ = -v^eV};} . For 
notational convenience, we remove the exponent 1 in the spherical tensors e 1 and r 1 ^. 
This should not bring confusion: if a vector takes part in a coupling, it is a first-rank 
tensor. Using the recoupling identity f|T3|) . we find the dipole transition intensity 

2 

\T fl \ 2 = 3{e*®r* fi }°{e®T fi }° = ^(-l)«{ e *Wr ■ {r}^r ft } a . (14) 

a=0 

Note that, for elliptic polarization, e is complex, whereas the wavevector k is always 
real. Each term of a decomposition over spherical tensors has often a clear physical 
meaning. In equation (fl~4"|) . the variables concerning the incident x-ray (i.e. e and e*) 
are gathered in {e*®e} a , the variables concerning the crystal are in {r^r fi} a . Thus, 
we can easily investigate the influence of a rotation R of the crystal on the absorption 
cross-section 

2 

a=0 
2 

a=0 

where we used equation ([2]) and that {r^^r/j} is a ath-rank spherical tensor. In 
particular, the spectrum of a powder sample is given by the average over all orientations, 
i.e. over all rotations R. This average is very simple with spherical tensors (D a (R)) = 
5 a)0 . Thus, the term a = gives the spectrum of a powder, called the isotropic spectrum. 

(IT^J 2 ) = {e* 8 e}° • {i>r /t }° = l(e* • e) (r* ft • Tfi ) = 

where we used equation f|T2l) and |e| 2 = 1. 

To interpret the term a = 1, we use the relation between vectors u, v and the 
corresponding first-rank tensors u 1 , v 1 (we restore the tensor rank in u 1 , v 1 for clarity): 
according to equation ( |27j) . {u 1 ®v 1 } 1 is the first-rank tensor corresponding to the vector 
(z/y / 2)u x v. Therefore, 

{e^e} 1 ■ {r}^} 1 = - ~(e* x e) • (r* fl x r fi ). 

The first cross product is related to the rate of circular polarization P c and to wavevector 
direction k of the incident x-ray by e* x e = —iP c k [5]. Moreover, the second cross 
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product is zero for a non-magnetic sample because time-reversal symmetry implies 
= Tfi. Therefore, the term a = 1 describes magnetic circular dichroism. 
The term a = 2 describes the linear dichroism of x-ray spectra. The number of 
non-zero components of {r^r^} 2 depends on the symmetry of the crystal [10]. We 
shall determine the number of non-zero components in the case of spinel and garnet. 



4-2. Quadrupole 

We consider the case of electric quadrupole transitions, which is more involved than 
dipole transitions. We start from the quadrupole transition operator T = e ■ r k • r and 
we rewrite it in terms of spherical tensors using equation (1281) : T = 3{{e®r} (g){k(g)r} } . 
In this expression e is coupled with r and k with r. As in the electric dipole, we want 
to gather all the terms concerning the crystal into a single tensor. For that purpose, we 
use equation ([TBI with the sum over g changed into a sum over a 

2 

T = ^(-l) a {e0k} a • {iw} a . 

a=0 

The term a = is zero because, according to equation (fl2|) . {e®k}° = — (1/ y/3)e ■ k = 
since the polarization and wavevectors are perpendicular. The term a = 1 is zero 
because equation ([27]) gives us {r^r} 1 = (i/y/2)r x r = 0. Thus, T is reduced to the 
single term 

T = {e®k} 2 • {iw} 2 = v / 5{{e®k} 2 ®{r®r} 2 } . 

The tensor {r®r} 2 can be expressed in terms of spherical harmonics (equation (23), p. 67 
of 0) 

where r, 9 and <fi are the spherical coordinates of r. For completeness, we give the 
components of {e®k} 2 : 

r ,-12 ( e z ^y) i^x i iky) 
{e®k} ±2 , 

2 _,_ ( e x i ity)k z + e z {k x i ik y ) 



{e®k} ±1 = =F 



{e®k} 



2 



x/6 V 2 

The electric quadrupole transition intensities are proportional to |7/j| 2 , where the 
transition amplitude is T fi = (f\T\i). Therefore, \T fi \ 2 = 5{P 2 ®Q 2 }°{R 2 ®S 2 } , with 
p2 = { e * 0k }2 ; Q 2 = (f\{ r ®r} 2 \i)*, R 2 = {e®k} 2 , S 2 = (f\{r®r} 2 \t) , 

4 

|T /t | 2 = 5{P 2 ®Q 2 }°{R 2 ®S 2 } = ^(-l) a {P 2 ®i? 2 } a • {Q 2 ®S 2 } a . (15) 

a=0 
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The apparent simplicity of this calculation is essentially due to the powerful tools given 
in reference [3]. A straightforward approach is quite heavy [10J. If the sample is 
nonmagnetic or the x-ray polarization is linear, the terms a = 1 and a = 3 are zero. 

As in the electric dipole case, the term a = corresponds to the isotropic spectrum 
obtained by measuring a powder. Equation (1291) gives us {P 2 ®i? 2 } = /c 2 /(2\/5) and 
the isotropic spectrum is 

The calculation of this average will be discussed in section 

If P 2 = R 2 (i.e. the x-rays are linearly polarized) or Q 2 = S 2 (i.e. the sample is 
non magnetic), then the terms a = 1 and a = 3 are zero. More generally, for any tensor 
T a with integer rank a, {T a ®T a } c is zero for c odd. This is due to the symmetry of the 
Clebsch-Gordan coefficients [21 E] (b(3aa\c~f) = (-l) a+b ~ c (aab(3\c-i): 

{T%T% = J2(^a(3\c 7 )T^ = {-l) 2a ^^ aaa ^ T ^ 
= (-1) C {TVT% 

where we first exchanged the summation variables a and (3, and used the symmetry 
of the Clebsch-Gordan coefficients, then the commutativity of T£ and TS, and the fact 
that a is an integer. 

If we consider only the case of linearly polarized x-rays or non-magnetic samples, 
only the terms a = 0, a = 2 and a = 4 are possibly non zero. The number of independent 
components depends on the crystal symmetry and is tabulated in [10] . 



5. Site symmetry 

We describe how to calculate the spherical tensor of a crystallographic site, assuming 
that it is invariant under the symmetries of the site. We can work in a reference frame 
corresponding either to the site symmetries or to the crystal symmetries. If we take the 
aluminum sites of spinels as an example, the site symmetry is whereas the crystal 
is cubic. In the reference frame adapted to the site symmetry, the z axis is along a 
diagonal of the cube. In the reference frame adapated to the crystal, the z axis is along 
an edge of the cube. The symmetrized tensors has usually less nonzero components in 
the site axes, but they are easier to calculate in the crystal axes. We also give the way to 
go from one reference frame to the other. Two examples are treated in detail: a spinel 
and a garnet. 



5.1. Symmetrized tensor 

If G is the symmetry group of the crystal, then a site has a symmetry group G' which 
is a subgroup of G. To know the general form of a spherical tensor invariant under the 
site symmetries, we start from a general tensor T and we calculate the symmetrized 
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tensor (T^) compatible with the site symmetries by using the classical formula 

= ]|m£ E T^jRi), (16) 

where i?' runs over the symmetry operations of the subgroup G' and where \G'\ is the 
number of elements of G' . From the physical point of view, this formula means that 
the tensor (T^) is obtained by averaging over all the symmetry operations that leave 
the site invariant. From the mathematical point of view, we project onto the subset of 
tensors that are invariant by any symmetry operation of G' . To check this, take any 
operation R in G' and evaluate the action of R on the symmetrized tensor: 

(RT^) = — — Tl l „D ( m i, m ,(R)D e m , m (R) 

R' m'm" 

= Tq>~\ l ^m" Dm"m {RR ) = (T^), 

' ' R' m" 

because, G' being a group, the set of operations RR' where R' runs over G' is the same 
as the set of operations of G'. 

5.2. Site and crystal axes 

The rotations R can be expressed either in the site axes or in the crystal axes. We shall 
see in the examples that the symmetric tensors are simpler in the site axes. Moreover, 
some computer programs need to be used in the site axesJJ However, the matrices R 
are easier to determine in the crystal axes, because they correspond to the symmetry 
operations of the crystal that leave the site invariant. Both cases will be treated in the 
examples of the following sections. 

It is also necessary to describe precisely how to go from one reference frame to 
the other. If ei, e%, are the orthonormal axes associated with the crystal and e' 1; 
e' 2 , e' 3 those of the site, we denote by R the rotation matrix such that e' { = RijGj. 
For example, if are the axes of the cube and are the trigonal axes along the 
( — 1,1,1) direction, then R is the rotation matrix corresponding to the Euler angles 
(0, arccos(l/\/3), 7r/4). It is the inverse of the rotation matrix of equation ( 1261) . It can 
be checked that e 3 = (— ei + e 2 + e 3 )/3 (i.e. the three-fold axis is the z-axis of the site 
along the (—1, 1, 1) direction of the cube) and e' 2 = (ei + e 2 )/ V2 (i.e. the y axis of the 
site is along the (1, 1, 0) direction of the cube). If R' is a symmetry operation in the 
cubic axes of the cube leaving the site invariant, the basis change formula in a vector 
space implies that RR'R~ l is a symmetry operation in the trigonal axes of the site. A 
spherical tensor T e will be denoted by T (3) when it is expressed in the trigonal axes 
and by T e (4) when it is expressed in the cubic axes. The argument 3 and 4 mean that 
the z axis is along a three-fold axis for a trigonal basis and a four-fold axis for a cubic 



| The most prominent example is the series of multiplet programs written by Cowan, Butler, Thole, 
Ogasawara and SearleplJ H2 UHl HH HS1 • 
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basis. The relation between T e (3) and T £ (4) is given by the formula 

7i(4)=^7i,(3)^ m (fl). (17) 

m' 

We give now two examples. 
5.3. The example of spinel 

We illustrate this method with the example of the aluminium site in spinel MgA^C^, 
which is the most common mineral of the spinel structural family. The spinel structure is 
derived from a face-centred-cubic close-packing of oxygens with a space group symmetry 
Fd3m. The conventional cubic cell contains 8 formula units, i.e. 32 oxygen atoms with 
24 cations in tetrahedral and pseudo-octahedral interstices. With origin choice 2 [16J, 
the Mg 2+ cations occupy 8 tetrahedral sites, which are located at the special 8a Wyckoff 
positions (1/8,1/8,1/8), with 43m (T^) point symmetry. The Al 3+ cations occupy 16 
pseudo-octahedral sites at the special 16d Wyckoff positions (0,1/4,3/4), with 3m (D 3( j) 
point symmetry. This symmetry corresponds to a small elongation of the octahedron 
along the trigonal axis, arising from a small departure of the position of the oxygens 
from the perfect fee arrangement. The O 2- ions are located at the Wyckoff positions 
32e (u,u,u) with point symmetry 3m. 



5.3.1. The site frame The simplest results are obtained when the reference frame of 
the site is used. The point group of the site we consider is D 3l j. The group D 3 d has six 
pure rotations and the same six rotations multiplied by the inversion. We assume that 
the property that we investigate is not sensitive to the inversion operator, so that we 
only have to consider the six pure rotations. It is natural to take the z-axis along the 
three-fold axis and the y axis along one of the Ci axes. The results do not depend on 
which C*2 axis is chosen. However, they would be different if the y axis were chosen, for 
example, between two Gi axes. The six pure rotations are: the unit, the C 3 rotation 
about the z-axis through the angle 27r/3, its square Cf, the Ci rotation about the y- 
axis through the angle n and the other two rotations C 3 C 2 and C|C 2 . These rotations 
have Euler angles (0, 0, 0), (0, 0, 2vr/3), (0, 0, 4vr/3), (0, vr, 0), (0, vr, 4tt/3) and (0, vr, 2vr/3), 
respectively. These rotations will be denoted by Ri, . . . , Rq. 

To calculate the symmetrized tensors for this site, we use equation (|T6|) . The special 
cases given in section 18.3.11 enable us to show that D^, m (i?i) = S m ' m , D t m , m (R2) = 
5 m > m e- 2mt7r/3 and D* m ,jR 3 ) = 5 m , m e~ im ™/\ Therefore, the sum £)J =1 D^jRj) is 
3^ m 'm if m is an integer multiple of 3, and zero otherwise. We calculate the Wigner 
matrices for the other three rotations and we obtain 

1 \ ^ n e ( t>i\ _ 1 \^ n e ( z> \ _ ^ m ' m + (~^Y m $m-m> 

1 1 R> j=l 

if m and m' are integer multiples of 3, and zero otherwise. Equation f|T6|) is then applied 
to a general fourth-rank tensor T^(3), where the argument (3) denotes the trigonal axes 
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and we obtain the non-zero components of the symmetrized tensor (T^(3)) 
Ctf(3)> = T 4 (3), 

<T 3 4 (3)> = - (T 4 3 (3)> = T34(3) /- 3(3) . 

For the second-rank tensor, all symmetrized components are zero, except for (T 2 (3)) = 
Tq(3). Of course, we have also the relation (Tq(3)) = T °(3), which is valid for any 
group. Now we show that time- reversal symmetry implies that the symmetrized tensors 
are real. According to equation (j3J) (To)* = Tq, so that Tq is real. Still by equation (J3]) 
we have (Tj)* = -T e _ 3 . Thus, (T 4 (3))* = (-T 4 3 (3) + T 4 (3))/2 = (T 3 4 (3)) is real as well. 

In x-ray absorption spectra, the symmetrized tensors are spectral functions 
depending on the photon energy. For the example of the electric quadrupole transitions 
we take, for each energy hu, 

T e = n 2 a Y,{(f\{^r} 2 \z)%(f\{mr} 2 \t)Y5(E f - E, - flu;), 
f 

where ao is the fine structure constant, Ei and Ef the energy of the initial and final 
states. The symmetrized tensors (T e (3)) can be calculated by multiplet programs. The 
value of these tensors for a chromium atom substituted for aluminum in spinel is given 
in figure [j] (see reference [1] for the details of the calculation) 



0.03 r 




-0.01 1 i i i i i i 

5986 5988 5990 5992 5994 5996 5998 

E(eV) 



Figure 1. The symmetrized tensors (T °(3)>, (T 2 (3)), (T 4 (3)) and (T 3 4 (3)) in the 
site frame, for the electric quadrupole transitions a the K-edge of a chromium 
atom substituted for aluminium in spinel. After averaging over the sites, only the 
symmetrized tensors (Tq ) and (Tq ) remain as independent parameters. 



5.3.2. The crystal frame We consider now the same average in the crystal axes. The 
chosen Al site with reduced coordinates (0,1/4,3/4) has a three-fold axis along the (- 
1,1,1) direction and a two-fold axis along the (0,1,0) direction. Therefore, the six pure 
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rotations of D 3d are now (i) the identity, denoted by (x,y,z), (ii) a C 3 rotation about 
(-1,1,1), denoted by (—y,z,—x), (iii) its square (—z,—x,y), (iv) a rotation of tt about 
(1,1,0) denoted by (y, x, — z), (v) a rotation of tt about (1,0,1) denoted by (z,—y,x), 
(vi) a rotation of tt about (0,1,-1) denoted by (— x, — z, — y). The notation used for 
the rotations is the result of the operation Rr in the cubic axes. For example, the C 3 
rotation gives 



Rr 



V 



1 
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\ 


/ 


X 


\ 




I 


-y\ 






y 








z 


/ 


V 


z 


J 




\ 


—X j 



The corresponding Euler angles are (0,0,0), (7r/2, ir/2, 0), (-7T, ir/2, 7r/2), (0,7r,7r/2), 
(0,7r/2, tt) and (3vr/2, vr/2, 3vr/2). 

We apply again equation (Tl6|) to the general second-rank tensor T^(4), where the 
argument (4) stands for the cubic axes. This gives us the symmetrized tensor (T 2 (4)) 

(r 2 (4)> =0, 

<T 2 2 (4)> = <T|(4)>* = -iX, 



with 



(^i(4)>= -<l?(4))* = (l + i)A, 
3f2f (4) - 3?T 1 2 (4) + ST 2 (4) 



A 



where we have used time- reversal symmetry as in equation ([3]). Note that A is real. For 
the tensor (T 4 (4)), 

= m (t 4 4 (4)) = y5 (T 4 4(4)) * = C; 



<T 4 3 (4)> = -(T 3 4 (4))* = (1-*)V7C, 
<r 4 2 (4)> = <T 4 (4)>* = 2zV2~C, 
(T\m= -(T 4 (4))* = (l + i)C, 



(19) 



with 



7T 4 (4) + V^O^ 4 (4) 



12 



C 



v/7(9fcT 3 4 (4) + 3T 4 (4)) + 2V23T 2 4 (4) + KT 4 (4) - 3T 4 (4) 



24 



Note that £ and £ are real. 



5.3.3. From site to crystal frame From this example, it is clear that the symmetrized 
tensor (T £ (3)) in the site frame is much simpler the same tensor (T e (4)) in the crystal 
frame. The relation between the trigonal and cubic axes in Butler's tables is worked 
out in section 18.41 of the Appendix. 
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To go from one to the other we apply equation (1171) and we obtain the relations 

= - mmi^mi, (20) 



{ T^ ))= i^2m^Mm. (21) 



We recover the relations (119p with 

2v/5(T 4 (3)) - Vl4(^(3)> 



c 



36 

7<T 4 (3)> + 2 V / 70(T 3 4 (3)) 



18 

For the second-rank tensor we find the relations (1TB]) with A = — (T 2 (3))/a/6. For 
the zeroth-rank tensor we have obviously (T° (4)) = (Tq (3)). 

5.^. T/ie example of garnet 

We consider now the Al site in garnet with the example of grossular Ca3Al2(SiOs)4, 
which is a cubic mineral, with the space group Ia3d. The cubic cell contains 96 oxygen, 
24 calcium, 24 silicium and 16 aluminium atoms. The Al 3+ cations are at the 16a 
Wyckoff positions (0,0,0). This site is a slightly distorted octahedral, with a small 
elongation along the (111) axis of the cube, and has the 3 (C 3i ) point symmetry. 

We calculate the symmetrized tensor in the site frame as for spinel, but with the 
smaller symmetry group C^. 

We find that the non-zero tensor components are [ID] (T 4 3 (3)), (Tq(3)), (T 3 4 (3)), 
(T 2 (3)) and (Tq(3)), as illustrated in figure [2] (see reference [17] for the details of the 
calculation). 

It would not be very illuminating to calculate directly the symmetrized tensor in 
the crystal frame. It is more interesting to rotate the tensor. Indeed, in the group 
C*3j, the z-axis is specified by the rotation axis, but the y-axis is arbitrary in the plane 
perpendicular to the rotation axis. This arbitrariness can be quite useful. For instance, 
the parametrization of the Hamiltonian is simplified by choosing the y axis so that a 
crystal field parameter is set to zero (see Ref. [12] p. 184). This simplifies the calculation 
of the eigenstates but the parameter reappears as the angle a between the y axis and 
the (1, 1, 0) direction of the cube in the plane perpendicular to the (—1, 1, 1) direction. 

The corresponding rotation matrix is 

^ cos(a + 7r/3) cos(a + 27r/3) cosa ' 

sin(a; + 7r/3) sin(a + 27r/3) sina 
V -1/V2 l/y/2 l/y/2 j 

The Euler angles are a, arccos(l/v^3), vr/4. The angle a describes a rotation about 
the axis (-1,1,1). Therefore, the a dependency of the result is very simple: we have 
Tm(a) = T^(0)e~ TOia because the corresponding Wigner matrix is D l m , m = S mm ie~ trna . 



R 2 
3 
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Figure 2. The symmetrized tensors (T °(3)), (T 2 (3)), (T 4 (3)) and (T 3 4 (3)) in the 
site frame, for the electric quadrupole transitions at the K-edge of a vanadium atom 
substituted for aluminium in grossular garnet. 



To calculate the symmetrized tensor (T^(4)) in the cubic axes, we use equation 
171) . we put s = (Tq (3)) and t r + iU = e~ ta (T^(3)) and we obtain 

70s + 20t r - Yli-fhi 
36 ' 
2\/35s - 7V2t r + 3iy/Eu 



{TtM)) = (I?(4)>* = 
<^ 3 (4)> = - <T 3 4 (4)>* = (1 - i) 
<T 4 2 (4)> = <T 2 4 (4)>* = i 



— 36 

TOs - V7t r 



(^(4))= -(T 1 4 (4))* = (l + ^) 



9 

2y/hs - y/Utr + 3iV42U 



36 

W( 4)> = - Zi±|^. (22, 

The symmetrized second-rank tensor in the cubic frame is the same as for spinel. 
The effect of the angle a on the experimental spectrum can be considerable, as is 
illustrated in figure [31 



6. From site symmetry to crystal symmetry 

We consider in this section another type of problem. We assume that we have calculated 
a symmetrized tensor (T £ ) for a certain site A. We want to know the value of the same 
tensor for all the sites equivalent to A. In the first section, we describe how this can be 
done. In an x-ray absorption measurement, we measure the sum of the signals coming 
from all sites of the crystal. Therefore, we do not really need to know the spectrum of 
each site, we need to know the spectrum averaged over all the sites. Two methods will 
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001 1 1 1 t i 1 j 1 i 1 1 1 

5470 5472 5474 5478 547B 5460 

Figure 3. The crystal symmetrized tensor (Tq ) as a function of a, for the electric 
quadrupole transitions at the K-edge of a vanadium atom substituted for aluminium 
in grossular garnet. 

be presented to do so: the coset method and the brute force method. Finally, we treat 
the examples of spinel and garnet. 

6.1. Changing site 

In general, the symmetries of the crystal make several equivalent sites. More precisely, 
if a crystal has the symmetry group G and the site has symmetry group G', then the 
number of equivalent sites in the crystal is the ratio |G|/|G"|. Assume that we have 
calculated a physical property described by a spherical tensor for a given site A. We 
want to calculate the properties of the equivalent site B. 

If ya and r# are the position vectors of site A and B, there is a symmetry operation 
R of G such that r B = Rya- If we denote by Ga and Gb the symmetry group of sites 
A and B, then the symmetry operation R' of Ga will be transformed into the operation 
RR'R' 1 of Gb- This can be checked because R transforms the rotation axis n of an 
operation R' of Ga into the axis -Rn, and RR'R -1 describes the rotation through the 
same angle as R' but about the axis Rn. This is indeed a symmetry operation of site 
B. Moreover, the tensor T^(B) at site B is related to the tensor T^(A) at site A by 
the relation 

1 

Tm{B) = T e m ,(A)D e rn , m (R x ). 

m'=-l 

We are now facing a typical subtelty of crystal symmetries. We could have expected 
the argument of the Wigner matrix to be R instead of -R -1 , but this is not the case 
because the rotation R is in fact a passive operation from the point of view of site B. 
The rotation transports the rotation axis, which is a vector of the reference frame. By 
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moving the atoms of the crystal, the rotation transports the reference frame of the site. 
Therefore, we are in the passive point of view and we need to use because our 
convention uses the active point of view. 

To check this, calculate the symmetrized tensor around B. 

(TL(B)) = r^- £ E T l'( B ) Dt m'm(R'B)- 
1 B| R' B eG B m'=-i 

Now, we can use the fact that the operations R' B of Gb can be obtained from the 
operations R' A of Ga by R' B = RR' A R~ l . Therefore 

(Ti(B)) = j±- E T UB)Dl, m {RR' A R-\ 

1 A| R' A aG A m'=-l 

where we used \G A \ = \Gb\, a consequence of the isomorphism between G A and Gb- 
We see that this is only compatible with the transformations T e (B) = T e (A)D e (R~ 1 ) 
and (T e (B)} = (T e (A))D i (R' 1 ). 



6.2. Changing the x-ray beam 

If we have calculated the spectrum of site A, it is possible to obtain the spectrum 
of any site equivalent to A by calculating the spectrum of site A for a rotated x-ray 
beam. This is physically clear because, if you rotate both the crystal and the x-ray, 
the spectrum does not change. Thus, if site A is measured with a polarization e and a 
wavevector k, the spectrum of the site obtained by rotating A with rotation matrix R 
is the same as the spectrum of A measured with a polarization i? _1 e and a wavevector 
i? _1 k. To be fully convinced, we prove that, if P a and Q a are two ath-rank tensors, 
then P a ■ (RQ a ) = {R^P a ) ■ Q a . Equations © and (PJ give us 

P a ■ (RQ a ) = P a ■ (Q a D a (R)) = Y,{-lVP-aQ a pDU R )- 

a/3 

The symmetry relation (see ref. [3], P- 79) D^ a (R) = (—l) 3 ^ a D a ^ a _ (3 (R~ 1 ) enables us to 
write 

P a ■ (RQ a ) = Ysi-lfP-aD-a-piR-^Q? 

a/3 

= ^(-l^p^ir 1 ))^ = (R^p*) ■ Q a . 

& 

For the example of the electric dipole transition probability (fbil) we find 

{e*®e} a ■ R({r* fl ®r ft } a ) = R-\{e*®e} a ) ■ {r* fl ®r fl } a . (23) 

It remains to prove that the rotation of {e*®e} a corresponds to the rotation of e* and 
e. This is done by using the following identity: 

R{P%Q b } c = {(RP a )®(RP b )} c . 



Site symmetry and crystal symmetry 



19 



To demonstrate the latter identity, we write the rotation in terms of the Wigner matrices 
and the coupled tensor in terms of the Clebsch-Gordan coefficients: 

R{P a ®Q% = Y,(™bP\ci)P:Q b pD^(R). (24) 

af/37' 

The classical identity ( !30l transforms this expression into 

R{P%Q b }c = (aa'bP'\ n )P^D a aaf (R)D b ^(R) 

a/3a'f3' 

= ^2(a^Wcry){RP)^(RQ)^ = {(RP a )r4RP b )Y r 

a'P' 

Equations ([23]) and fl23D yield 

{e*®e} a • R({r* fi ®r fi } a ) = {(R- 1 e*)®(R- 1 e)} a ■ {r* fi ®r fi } a . 

In other words, rotating the crystal by R gives the same result as rotating the x-ray 
beam by R" 1 . The same result is true for the electric quadrupole transition probabilities, 
except for the fact that the polarization k are rotated by R~ x and not only e. 

6.3. The coset method 

The coset method is a powerful way to calculate the the tensor symmetrized over 
the crystal from the tensor symmetrized over a single site. We first introduce some 
mathematical concepts. 

6.3.1. Mathematical aspects The symmetry group Ga of a site A is a subgroup of the 
symmetry group G of the crystal. The number of sites equivalent to A in the unit cell 
is the ratio (GI/IG^I, where \G\ and \Ga\ denote the number of elements of G and Ga- 
This ratio, denoted by n in the following, is an integer by the Euler-Lagrange theorem 
([7j p. 11). We can pick up n symmetry operations gi in G that transport site A to the 
n sites equivalent to A. 

We first introduce some mathematical concepts [7] . If g is an element of a group G 
and Ga a subgroup of G, the set gGA = {gh : h G Ga} is called a coset. If we take two 
elements g and g' of G, then either gGA and g'GA are either identical or disjoint (i.e. 
they have no element in common). The number of different cosets is n = \G\/\Ga\ and 
the number of elements in each coset is \Ga\- Moreover, every element of G belongs to 
one and only one coset. Therefore, if we pick up any element $ in each coset, we have 
G = g\GA U . . . U g n GA and each g^ is called a representative of its coset. 

6.3.2. Cosets in a crystal We apply these concepts to a crystal. The crystal symmetry 
group is a space group G. The symmetry group of a site A is the set of operations of G 
that leave site A invariant. More precisely, if is the coordinate vector of site A, then 
Ga = {g £ G\g(rA) = r^}. It is clear that Ga is a subgroup of G. If n = \G\/\Ga\, 
there are n sites in the crystal that are equivalent to A. The group G is partitioned 
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into n cosets g\G^ . . . , g n GA- All the elements of a given coset send site A to the same 
equivalent site. Take an element #j in each coset giGA- 

The symmetrized tensor (T e )x over the full crystal is obtained from the site- 
symmetrized tensor (T )a by the operation 

n 

( Tt )x = -Y / ( T 'U Dt i9T 1 ), (25) 
i=i 

where, for a space group operation g, D e (g) is the Wigner matrix corresponding to the 
proper rotation of g. Equation (I25I) can be described by saying that we average over 
the symmetries of site A, then we go to another sites with g^ and we average over the 
symmetries of that other site. 

We show now that equation (1251) gives the same result as an average over all crystal 
symmety operations. We already know that 

(T% = £ T*D%) = -L T i D e (h~ 1 ), 

' A ' heG A ' A ' HeGa 

because the sum over the elements of a group is the same as the sum over the inverse 
elements of this group. Therefore 



1 n 1 n 

' A ' i=l h&G A ' ' i=l heG A 



IGI ^ ^ ' ' \G\ 

1 1 i=i heGA g&G 

where we used (gih)^ 1 = h^g^ 1 and |G| = n|GU|- Note that the proof holds because 
we used sets g«GU (i.e. left cosets) and not G^ft (i.e. right cosets). 

In other words, the average over the site symmetries followed by the average over 
the sites gives the average over the crystal symmetries. This can be considered as 
a factorization of the average, because ^ g D l (g) = '52 h D e (h)^2 i D e (g^' 1 ). It can be 
checked easily that the result of equation (}2"5]) is the same if we replace gi by any g\ in 
giGA- We illustrate this by our two favorite examples, spinel and garnet. 

6.3.3. The example of spinel We call site A the site (0,1/4,3/4) of the spinel structure. 
The space group G has 192 operations^. Twelve of them leave site A invariant: 
(x, y, z), {-y + 3/4, z, -x + 3/4), {-z + 3/4, -x + 3/4, y),(y + 1/4, x + 3/4, -z + 1/2), 
(z + 1/4, -y + 1/2, x + 3/4), (-x, -z + 1/2, -y + 1/2), (-x, -y + 1/2, -z + 1/2), 
(y + 1/4, -z + 1/2, x + 3/4), (z + 1/4, x + 3/4, -y + 1/2), {-y + 3/4, -x + 3/4, z), 
(—2 + 3/4, y, — rr + 3/4), and (x, z, y). This set of twelve operations is a group isomorphic 

§ The full space group has of course an infinite number of operations because it contains all the 
translations by a vector of the Bravais lattice. We consider here the group G obtained as the quotient 
of the full space group by the group of translations of the simple cubic lattice. We could have used 
the quotient of the full space group by the group of translations of the fee lattice (which has only 48 
operations) , however our choice is simpler from the programming point of view: we take the symmetry 
operations given by the table [16) . and we identify two points whose coordinates differ by an integer. 
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to D 3 d- The isomorphism ip is described explicitly as follows. If is the coordinate 
vector of site A, for any operation g of the set, we define the operation <p(g) by 
<p(g)(r) = g(r + r^) — r A . It is an isomorphism because <p(g'g) = (p(g')ip(g): 

^<yM#)( r ) = ^<y)(#( r + *a) - r A ) = (g'g)(r + r A ) - r A = ip(g'g)(r), 

and it can be checked that the image by ip of the first six operations are the rotations 
listed in section 15.3.21 the image of the other six operations are the same rotations 
multiplied by -1. The space group operations g are selected by the condition that the site 
is fixed: g(v A ) = r A . Thus, the origin is a fixed point of <p(g): ip (g)(0) = g(r A ) — r A = 0. 

There are 16 sites equivalent to A because the ratio \Fd3m\/\D 3d \ is 16. However, 
each site is equivalent to three other sites by pure lattice translations (1/2,1/2,0), 
(1/2,0, 1/2) and (0, 1/2, 1/2). The x-ray spectrum of these sites will be equal because 
their orientations with respect to the x-ray beam are the same. Therefore, we are left 
with 4 equivalent sites: A itself and the sites with coordinates (1/4, 3/4, 0), (3/4, 0, 1/4) 
and (1/2, 1/2, 1/2). A representative of the coset corresponding to each of these sites is 
(-y + 1/2, x + 3/4, z + 1/4), (-x + 3/4, -y + 3/4, z) and (y + 1/4, -x, z + 1/4). 

6.3.4- The example of garnet We consider now site A as the Al site in grossular garnet 
with the three-fold axis along the (—1, 1, 1) direction. Its coordinates are (1/2,1/2,0). 
It is left invariant by the following six operations: identity: (x,y,z), rotation through 
27r/3 about the (—1, 1, 1) axis: (—y,z + 1/2, — x + 1/2), rotation through 47r/3 about 
the (—1, 1, 1) axis: (—z + 1/2, — x, y + 1/2), and the same operations multiplied by an 
inversion (—x, —y, —z), (y, —z + l/2,x + 1/2) and (z + 1/2, x, —y + 1/2). This group is 
C*3j. The number of cosets (i.e. the number of sites equivalent to A) is |/a3rf|/|C3i| = 16. 
If we remove the lattice translation (1/2,1/2,1/2) we are left with 8 equivalent sites. 

The approach in terms of cosets is quite powerful in practice because it completely 
avoids the explicit description of the sites and of the symmetry operations that transform 
a site into another one. Many programs compute the symmetry operations of the sites 
(for example PWSCF). It is then enough to take one of these sites, to choose any 
representative gi in each coset and to calculate the contribution of all equivalent sites by 
the formula (T e ) A D £ (g^ 1 ). The average over the crystal is then obtained with equation 

(ESD. 

6-4- The brute force method 

If one is not interested in the contribution of each site to the spectrum of the crystal, a 
still simpler solution is to take the average of (T e ) A over all the symmetry operations of 
the crystal. This is not very clever because the site operations have already been taken 
into account and we average over them a second time, but this method is well suited 
for computers. 

We show now that averaging the site-symmetrized tensor over all the symmetry 
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operations of the crystal gives the same result as the coset method. 

±J2(T%D\ 9 ) = -\- J2 J2T*D\h)D\g) 

' ' geG \ A\\ | heGA g£G 

= ^7p E T,T'D'(k)D'(h')D'( 9 ^ 

1 1 h,h'<=G A i 



1 711 heG A i 

where we used the identity J2hh>eG A -D e (hh') = \Ga\ J2heG A D e (h), which is readily 
established. Therefore, the average over all the symmetry operations of the crystal gives 
the same result as the average over the sites, irrespective of the number of equivalent 
sites. 

6. 5. The case of spinel 

We illustrate the coset method with the case of spinel. We first notice that the absorption 
cross section is invariant under a translation of the Bravais lattice, because such a 
translation multiplies the wavefunction by a phase (independent of r) that disappears 
in the square modulus. Therefore, by removing the translations, we can replace the 
representatives of the four cosets given in section 16.3.31 by the four rotations about the 
z-axis of the crystal through angles 0, 7r/2, 7r and 3n/2. For a fourth-rank tensor, the 
average over the coset operations is rather drastic. The only non-zero elements of the 
matrix M = (1/4) £\ D A (gr 1 ) are M(-4, -4) = M(0, 0) = M(4, 4) = 1. Therefore, the 
crystal-symmetrized fourth-rank tensor is (T A )x = (T 4 (4))M: 

1 V 4 ) x = (T 4 (4)>. 

The relation between (Tj 4 (4)) and the site-symmetrized tensor in the trigonal axes is 
given by equation fT20l) . 

For a second-rank tensor, the matrix M = (1/4) ^2i i D 2 {g~ 1 ) has a single non-zero 
element: M(0, 0) = 1. Therefore (Tq)x = 0, as expected [10J. 

6. 6. The case of garnet 

We illustrate the brute force method with the case of garnet. To calculate the spherical 
tensor of garnet, we use the brute force method and calculate M = (1/48) J2rD £ (R), 
where the sum runs over all the symmetry operations of the cube. 
For I = 4, the only nonzero matrix elements are 

M(-4, -4) = M(-4, 4) = M(4, -4) = M(4, 4) = A 




70 

M(-4, 0) = M(0, -4) = M(0, 4) = M(4, 0) = ^ 
M(0,0) =1 
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Therefore, the only non-zero components of a fourth-rank tensor are 



(T*} x = \-{Tt) x = \hr(T\)x = (T 4 (4)> 



5 V 5 ' vy " 18 

where we recall that s = (T 4 (3)) and t r = (l/2)( e - 3ia (T 3 4 (3)) - e 3ia (Ti 3 (3))). Note that 
s and t r are real because of time-reversal symmetry. 

7. Spherical average 

In Cartesian coordinates, the electric quadrupole matrix elements can be written, for 
linearly polarized x-rays 

ijlm 

with (Tijim = (i\xiXj\f)(f\xix m \i). For a powder sample we have the spherical average 

) ~\~ ^^"yyyy 2o~ zzzz &0~xyxy &0~xzxz ^O'yzyz 

0~xxyy ®~xxzz @yyxx 0~ zzxx 0~yyzz 0~zzyy)- 

If the system is non magnetic, then o~ajj = ajja and the average further simplifies 

) (^xxxx ®~yyyy "zzzz 3(J X yxy 3o~ X zxz 3(Jy Z y Z 

0~xxyy ®~xxzz @yyzz)- 

8. Formulas 

8.1. Rotation matrix 

The rotation of ip around the direction n (a unit vector) is represented by the rotation 
matrix (see Ref. [2], p. 10) R — Id + sinipN + (1 — cos?p)N 2 , where N is the skew- 
symmetric matrix with matrix elements = — e^nk, so that (N 2 )^ = UiUj — Sij. 

Conversely, the rotation angle ip and the rotation axis n are determined from the 
rotation matrix R by the relations cosip = (tr R — l)/2, nisini/j = (R32 — -R23V2, 
r^sin^ = (#13 — i? 3 i)/2, and nssinip = (R 2 \ — R\2)/2. This is a corrected version of 
the relation given in Ref. [2] p. 20. 

8.1.1. Euler angles Ref. 0, p. 24 

Cqi dy S Q^Siy S-y S C-y Cqi S ^ 

R = 

where < a < 2tt, < (3 < tt and < 7 < 2n, c a = cos a, s a = sin a, etc. There is a 
one-to-one correspondence between rotations and parameters in this range, except for 
the case (3 = and (3 = it, which describe the rotation through the angle a + 7 and 
a — 7, respectively, about the axis (0, 0, 1). 
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8.1.2. Euler-Rodrigues parameters Although the Euler angles are more common, the 
Euler-Rodrigues have the advantage of being expressed in terms of the rotation matrix 
without using circular functions. Thus, it can be convenient to derive analytical 
expressions. From the axis n and the angle if), we define the Euler-Rodrigues parameters 
(Ref. [2J p. 54). a = cos(if>/2), on = sm(if>/2)rii. These parameters are interesting for 
analytic calculations because they can be derived from the rotation matrix without 
using trigonometric functions ([2] p. 19): If tr R ^ — 1, then ao = y/tr R + 1/2, 
"i = (#32 - J R 2 3)/(4a ), a 2 = (R13 - i? 3 i )/(4ao), an d a 3 = (R 21 - R 12 )/(4a ). If 
tr R = —1, then a = and = (signa i )y / (l + Ru)/2 for i = 1, 2, 3 with signer = 1, 
signer = sign_Ri2, signa 3 = signi?i 3 . 



' «q + a\ — a 2 — a\ 2ct\a 2 — 2ao«3 2a\a 3 + 2otQOt 2 

R = 2a±a 2 + 2aoa 3 «q + a 2 — a\ — a\ 2a 2 a 3 — 2a^ai 

i 2ai« 3 — 2a a 2 2a 2 a 3 + 2a «i cyq + oi\ — a\ — af i 

where a = cos(^/2), a, = sin(^>/2)nj. 



8.2. Solid harmonics 

For a vector r = (x,y, z), the solid harmonics Y^iv) are defined by (see Ref. [2] p. 71) 

\2l + 1)(£ + m)\(£ - to)! (~x - %) fe+m (x + i y ) k z £ - 2k - m 
4^ 2 2fc+m (fc + m)\k\(l -m- 2k)\ ' 



Y e m r 



where runs from max(0, —to) to the integer part of (£ — m)/2. The most important 
example of solid harmonics is 




3 

47T 



/ 



x—iy 
Z 



\ 



i _x±iy . 

\ y/2 J 



where the upper component is Y 1 (r). 



5.5. Wigner matrices 

There are several representations of the Wigner rotation matrices. We present here the 
expressions in terms of Euler angles and of Euler-Rodrigues parameters. Other formulas 
have been derived, for example the recent invariant spinor representation [T8] . 

5.5.1. Euler angles For a rotation R expressed in terms of Euler angles a,/3,j, the 
Wigner matrix is ([2] p. 46) 

Dmi m {R) = e %m a d^ n , m (f3)e 11711 . 

Various expressions exist for the reduced Wigner matrix d^, m (/3). The following formula 
(valid for half-integer £) is particularly convenient for computers, because it avoids the 
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presence of singular terms (|2| p. 50): 



where = min(£ + m, £ — m, £ + m', £ — m') and the non-negative integers /i, z/ and A 
are determined by the value of k 

• if k = £ + m, then fi = m' — m, u = —m! — m and A = m' — m', 

• if = £ — m, then fi = m — m', v = m! + m and A = 0, 

• if k = £ + m', then /i = m — m', z/ = —ml — m and A = 0, 

• if k = £ — m', then /i = m! — m, v = —w! + m and A = m! — m'. 

In this expression, the only possible numerical difficulty occurs for 0° that should be 1. 
The Jacobi polynomials Pj: (x) are given by the formula 

^(^tCr)(r_i)(^)"(T)' 

For example, the Wigner matrix for first-rank tensors is 



D 



/ cos/3+1 i( a +y) sin/3 iq cos/3-1 j( a - y ) \ 

' 2 ^ y/2 ^ 2 ^ \ 

-^e ! ^ cos/3 ^e-*T 

. cos/3-1 j(7-q) sin/3 -jq cos g+l -t(q+ 7 ) , 

V 2 ^ C 2 ^ / 



where the upper left matrix element is D\ x _ x . Two useful special cases are d^ m ,(0) = 
cW and ^ TO /(7r) = (-\Y +m 5 m - m >. 

8.3.2. Euler-Rodrigues parameters The Wigner rotation matrix is 

D L>m( R ) = + m')\(£ - m')\{i + m)\(£ - m)\ 

^ (a - ia^Y +m - h {-iai - a 2 ) m '- m+k (-ia 1 + a 2 ) k (a + ia 3 )^ m '- fc 
^ (£ + m - k)\(m' - m + k)\k\(£ - rri - k)\ ' 

where k runs from max(0, m — ml) to min(£ + m, £ — m'). 
8.4- Butler's orientation 

The powerful multiplet program developed by Theo Thole and Barry Searle is based on 
Butler's conventions. For the calculation of trigonal sites in cubic crystals, it is necessary 
to know precisely the relation between the cubic and trigonal reference frames, which is 
not clearly stated in Butler's book. To determine it, we combine Butler's tables pp. 522, 
527 and 549 of [12J. This shows that the transition between spherical harmonics |lm) 3 
in the trigonal axes (i.e. in the O-D3-C3 basis) and spherical harmonics |lm)4 in the 
cubic axes (i.e. in the O--D4-C4 basis) is 

H - i) 3 = n - ^(1-0(^+1) + , 10> 1 + |n)i (i + 0(VS - 1) 



24 V24 
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|10>; 



1)^ + 110)^ + 111). 



1 + i 
V6 ' 



i) 4 (1 -^- 1) + |io) 4 ^ + |ii) 4 ^^ +1) 

v 24 v 3 
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|H>3 

This can be rewritten 

i 

m'=—l 

for the rotation R corresponding to the Euler angles a = Sir/A, j3 = arccos(l/v^3) and 
7 = 7r. This corresponds to the C 3z axis of D 3 along the (-1,1,1) direction of the cube 
and the Ci y axis of D 3 along the (1,1,0) direction of the cube (see fig. 11.6 of Ref. [12J, 
p. 204). The inverse rotation has Euler angles (0, j3 , vr/4). To be still more detailed, 
the rotation 

/ l/y/6 1/V2 -l/y/3 \ 
R = i?(37r/4,/3 ,7r) = -1/VQ l/y/2 l/y/3 , (26) 

V 1/V3 / 

transforms any symmetry operation R' in the D 3 d axes into the symmetry operation 
RR' R~ Y in the cubic axes. 



9. Coupling identities 

We gather some useful coupling formulas. If a, b, c and d are vectors, we denote by 
a 1 , b 1 , c 1 and d 1 the corresponding first-rank spherical tensors. Then, according to 
reference [3] p. 66 and 67, 

{aW} 1 = -5=(a x b) 1 . (27) 
v2 

{{aW} Vc^d 1 } } = i(a ■ b)(c ■ d). (28) 

{{aW} V^d 1 } 2 } = (a ■ c)(b ■ d) + ~(a • d)(b ■ c) 

-i(a-b)(c-d)). (29) 
To prove equation f[T3"j) . we start from the identity 

ah c 



{p a ®Q h Y ■ {R d 9 s e y = (-i) 2a+b - d J2( 2c 



e d g 



9 

x {P a ®R d y ■ {Q b ®S e } 9 , 

where g runs between max(|a — d\, \b — e|) and min(a + d, b + e) (see equation (13) p. 70 
of [3]). Equation (TlBl corresponds to the case c = because of the special value of the 
6j-symbol (equation (1) p. 299 of [3]) 

a b I _ (_l)a+e+f $abO~de 

edgj ' v /(2a + l)(2d+l)' 
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The interplay between Wigner matrices and Clebsch-Gordan coefficients is 
described by the following identity (equation (5) p. 85 of [3]) 

J2(aab(3\ C1 ')D c ^(R) = ^ (aa'b(3' \ ay ) D a aa , (R) Dpp, (R) . (30) 

7' a' /3' 
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